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We calculate the shear viscosity ri ~ rjcfi+rin in a neutron star core composed of nucleons, electrons 
and muons (770^ being the electron-muon viscosity, mediated by collisions of electrons and muons 
with charged particles, and rjn the neutron viscosity, mediated by neutron-neutron and neutron- 
proton collisions). Deriving rje^, we take into account the Landau damping in collisions of electrons 
and muons with charged particles via the exchange of transverse plasmons. It lowers »7e^ and leads 
to the non-standard temperature behavior r/eti oc T~^^^ . The viscosity r^n is calculated taking into 
account that in-medium effects modify nucleon effective masses in dense matter. Both viscosities, 
00 , Vcfj. and rin, can be important, and both are calculated including the effects of proton superfluidity. 

. They are presented in the form valid for any equation of state of nucleon dense matter. We analyze 

' the density and temperature dependence of 77 for different equations of state in neutron star cores, 

I and compare 77 with the bulk viscosity in the core and with the shear viscosity in the crust. 

PACS numbers: 97.60.Jd, 52.25.Fi, 52.27.Ny 
,—1 ■ I. INTRODUCTION 

<N : 

^ Neutron stars are very compact. Their typical masses are ^ IAMq (where Mq is the mass of the Sun), while their 

, radii are as smaU as 10 km. As a result, a neutron star core contains matter, whose density p reaches several po 
(po ~ 2.8 X 10^"* g cm~'^ being the density of the standard saturated nuclear matter). The core is composed of uniform 
neutron-rich nuclear matter and extends from p w 0.5 po to the stellar center (where p can be as high as lOpo)- It 
attracts special attention because of its poorly known composition and equation of state (EOS); e.g., Ref. From 
' outside, the core is surrounded by a thin (~ 1 km thick) and light (a few per cent by mass) crust composed of atomic 
d . nuclei, strongly degenerate electrons and (after the neutron drip at p > 4 x 10"'^^ g cm~'^) free neutrons. 

In this paper, we study the shear viscosity of neutron star cores. It is an important transport property which 
, affects the relaxation of hydrodynamic motions, particularly, a possible differential rotation within the star and stellar 
^ ' oscillations The shear viscosity can be important for damping gravitational wave driven instabilities (for instance, 
00 \ r- modes; e.g., and references therein). Its knowledge is required to analyze the efficiency of such instabilities for 
' generating gravitational waves. 

' For simplicity, we consider the cores composed of strongly degenerate neutrons (n) , protons (p) , electrons (e) , and 
. muons (/i) - npe/i-matter, neglecting a possible appearance of hyperons and/or exotic forms of matter (pion or kaon 
00 condensates or quarks or their mixtures) as predicted by some EOSs at p > 2po; see, e.g., Ref. jlj]. The electrons 
, and muons constitute almost ideal gases. The muons are absent in the outermost part of the core. They appear at 
• densities exceeding a threshold value P/i ~ po 01 at which the electron chemical potential reaches the muon rest-mass 
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energy (/Xc = rrif^c w 207 rricc ). The electrons are ultra-relativistic, while the muons are non-relativistic just after 
the threshold but become relativistic at higher p. In contrast to electrons and muons, nucleons constitute a strongly 
interacting Fermi liquid where protons are essentially non-relativistic, while neutrons become mildly relativistic at 
p > 2po. The neutrons and protons can be in superfluid state (e.g., Ref. (H). 
■ The main contribution to the shear viscosity 77 in a neutron star core comes from electrons and muons (lightest and 
most mobile particles) and neutrons (most abundant particles), 

■q^Vci^ + Vry- (1) 

The viscosity 770^ of electrons and muons is mainly limited by collisions of electrons and muons between themselves and 
with other charged particles (protons, in our case) via electromagnetic forces. In contrast, the neutron contribution r/n 
is limited by neutron-neutron and neutron-proton collisions mediated by strong interactions. As a result, 770^ and ?7n 
are nearly independent (belong to different - electromagnetic and nuclear - sectors) and can be calculated separately 

In applications, one often employs the viscosity 770^ calculated by Flowers and Itoh [6] for non-superfluid matter. 
Recently, Andersson et al. j?! have estimated rj^^ for superfluid matter. However, these studies neglect an enhancement 
of collisions of relativistic charged particles due to the exchange of transverse plasmons. The significance of this effect 
was demonstrated by Heiselberg and Pethick l8| in their study of transport properties of ultra-relativistic quark 
matter. Recently we (Shternin and Yakovlev [9| - hereafter SY07) have reconsidered the electron-muon thermal 
conductivity taking into account the exchange of transverse plasmons. This effect can reduce by several 
orders of magnitude. 
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Here we reanalyze rjofi in the same manner. We closely follow SY07 and omit technical details. In addition, we 
reconsider ?7„, which is a more difficult task involving nucleon-nucleon collisions. The viscosity rj^ was calculated by 
Flowers and Itoh Q for one EOS of non-superfluid matter assuming in-vacuum nucleon-nucleon scattering. These 
results were fitted by Cutler and Lindblom 2] by a simple analytical expression which is widely used; according to 
Ref. @, 7?n > Vcfj.- Recently Bcnhar and Valli 10] have calculated yyn for pure neutron matter in a self-consistent 
manner using the same nucleon interaction potential to derive rjn and construct the EOS (also for one EOS). We 
calculate ?7n in a more general way than Flowers and Itoh [6^. Our approach is similar to that used by Baiko, Haensel 
and Yakovlev [llj (hereafter BHYOl) for evaluating the thermal conductivity of neutrons. In addition, we employ 
recent developments [l^] in calculations of nucleon-nucleon scattering cross sections in nuclear matter. As in BHYOl, 
we take into account superfiuidity of protons. Again, we closely follow the derivation of BHYOl and omit the details. 

After calculating r/e^ and ?7n, we analyze the shear viscosity in neutron star cores with different EOSs. 



II. SHEAR VISCOSITY IN NON-SUPERFLUID MATTER 



The shear viscosity is calculated from a system of coupled Boltzmann kinetic equations 

i 

where Fc is the distribution function of momentum-transfer carriers c (with c = e, yit, or n, in our case); i =n, p, e, fj, 
runs over all particle species; Vc is the velocity of particles c, and Id is a collision integral, that describes a scattering 
of particles c and i: 



I = I V 



J dp2dpydp2' W„(12|1'2') 



X [Fi,i^2'(l - i^i)(l - F2) - FiF2il - FiO(l - i^2')] ■ (3) 

Here, 1 and 2 denote particle states before a collision; 1' and 2' are particle states after the collision; p is the particle 
momentum, a is the spin state, and Wd is the differential transition probability. The Kronecker delta Sd is included 
to avoid double counting of collisions between identical particles {c = i). 

Distributions Fc slightly deviate from the equilibrium Fermi-Dirac distributions fc owing to the presence of a small 
hydrodynamical velocity field Y, 

Fc = fc-'^.^J^, fc^Uvi'-^]+l\\ (4) 



dec' ' \ \ keT , 

where Sc is the particle energy, ^Lc is its chemical potential, T is the temperature, ks is the Boltzmann constant, 
and measures a deviation from equilibrium. The electron-muon and neutron transports are decoupled because 
we neglect electromagnetic interaction between the leptons and neutrons. For calculating "q^^^ the electrons and 
muons are treated as the only momentum carriers which undergo collisions between themselves and with protons. For 
calculating 77,1, the only momentum carriers are assumed to be neutrons, while the contribution of protons is neglected 
due to their small fraction. Therefore, the protons are thought to be passive scatterers which obey the equilibrium 
Fermi-Dirac distribution. Nonequilibrium parts of the electron, muon, and neutron distributions are found using the 
standard variational approach with the simplest trial function, 

'^c = -Tc [vcaPcP ~ VcPc5al3^ Vaft, (5) 

where Tc is an effective relaxation time of particles c, Vc is their velocity, and 

with ■^"^ = divV = 0. 

The resulting shear viscosity is expressed through the effective relaxation times in a standard way, 

ncPp^Tc IT'f.Pl'f.Tf, UnPp^Tn , , 

om* 5m* 5m* 
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where rje, rj^j^, and rj^ are, respectively, the partial electron, muon, and neutron shear viscosities; ric is the number 
density of particles c; ppc is their Fermi momentum; and m* is an effective mass on their Fermi surface. The electron 
and muon effective masses differ from their rest masses due to relativistic effects, m* = fic/c^ and m* = A'p/c^- The 
neutron and proton effective masses differ from their bare masses mainly due to many-body effects in dense matter 
(being determined by neutron and proton densities of state near appropriate Fermi surfaces). 

Linearizing the kinetic equations, multiplying them by [viaPip — ■|wiPi<5q/3), summing over ci and integrating over 
(27r?i)^'^dpi, we obtain a system of equations for the relaxation times. 



e, fi, n, 



(8) 



where we introduce the effective collision frequencies, 

in'^h^ f dpi dpi> dp2 dp2' 



2 4 I -'^ 2 2 



(27rft)i2 

-(Pi-Pi')' 
dpi dpi' dp2 dp2' 



(27rft)i 



l^,,(12|l'2')/i/2(l-/i')(l-/2' 



M^„(12|l'2')/i/2(l-/i')(l-/2') 



^22 I22,/ \2 / \2 

3P1P2' - 3P1P2 + (Pl •P2) -(Pl-P2') 

with tyci(12|l'2') = (1 + (5cj)"^ J2 Wci(12|l'2') (the sum is over spin states of ah particles 1,2,1', 2'). 

spins 

The formal solution of ^ for the npe^-matter is 



(9) 



(10) 



where 



C/J, 



/J.C 



(11) 



In the absence of muons, the expression for rjofi simplifies, 



(12) 



Once collision frequencies are found, the viscosity is obtained from Eq. ([7]). In order to determine the collision 
frequencies from Eqs. ^ and PU]) one needs to know the transition probability Wci(12|l'2'). The colhsions of 
charged particles should be considered with a proper treatment of plasma screening of electromagnetic interaction. 
We discuss the plasma screening and the calculation of 77,, in Sees. Ill AHII Cl The neutron viscosity Tyn is studied in 
Sec. IIIDl In Sec. [Illwe consider nonsuperfluid nucleons; the effects of proton superfluidity are analyzed in Sec. IIIII 
Throughout the paper we use the simplest variational approach. A comparison with an exact solution is discussed in 
Sec HTUH 



A. Plasma screening 

The plasma screening in neutron star cores is discussed in SY07. Here, we outline the main points. 
The differential collision probability can be written as 

T^„(12|1'2') = 4^^^ 5{p, +p2~ pv ~ P2')Siei + S2 - er - £2') t^^, (14) 

il 1 -\- Oci 
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where {\Mci\ ) is the squared matrix element summed over final and averaged over initial spin states. For collisions of 
identical particles, we have Mcc = Mc^^ —Mcc\ where the first and second terms correspond to the scattering channels 
12 —^ 1'2' and 12 2'1', respectively. Collisions of different particles go through a single channel, Md = M^l\ 

where hq — — pi and fku — ey — ei are, respectively, momentum and energy transfers in a collision event; 

"^c'c ~ ("^c^l y "^c c) = (2TO*c)~^(uc'7'^itc) is the transition 4-current {v = 0, 1, 2, 3), Jtc'c is the component of Jc'c 
transverse to q; j'^ is a Dirac matrix; Uc a normalized bispinor (with UcUc = 2racC^), and Uc is a Dirac conjugate. 
The first term in Eq. (|15p corresponds to direct Coulomb interaction via the longitudinal currents (with respect 

to q) ; the space-like longitudinal component of the current is expressed through the time-like component J^^l^ with 
the aid of charge conservation condition. The second term describes the interaction via transverse currents. It is 
especially important for relativistic particles because Jc'c/J^'^i ^ Pc/ {fTi*cC)- Longitudinal and transverse interactions 
are accompanied by different plasma screenings described by the functions lit and 11; in the denominators of Eq. (jlSp . 

The collision energy and momentum transfers in neutron star cores are typically small, huj ^ ksT <C Si and 
hq <^ pFi. This smallness allows us to use the weak-screening approximation which greatly simplifies the consideration. 
Moreover, one typically has ui <C qvpi, so that it is sufficient to use the asymptotic expressions (e.g., SY07) 

2 4a ■r-^ 

n; = qi = Z^m^PFiC, (16) 

i 

4 qc fi^ qc 

where a = e^/hc ~ 1/137 is the fine structure constant; qi and qt are characteristic plasma wavenumbers which 
depend on plasma composition (summation is over all types of charged particles); qi is the familiar Thomas-Fermi 
screening wavenumber; qt < qi, with qt qi in the limit of ultra-relativistic particles. Longitudinal interactions 
(via the exchange of longitudinal plasmons) are mediated by static non-dissipative screening with characteristic 
wavenumber qi (11; is real) , while transverse interactions (via the exchange of transverse plasmons) are accompanied 
by the coUisionless Landau damping (IIj is purely imaginary). Characteristic momentum transfers in transverse 
interactions arc A = (7rLj/(4:cqt)Y^^qt qi, meaning that such interactions occur on larger spatial scales than 
the longitudinal ones. Therefore, for relativistic particles, the transverse interactions can be more efficient. The 
importance of such interactions was pointed out by Heiselberg and Pethick Q in their study of kinetic properties of 
relativistic quark plasma. So far in all calculations of kinetic properties in neutron star cores (except for SY07) the 
transverse interactions have been erroneously screened by the same static dielectric function 11; as the longitudinal 
interactions. This approximation strongly (up to several orders of magnitude) overestimates the electron- muon thermal 
conductivity (SY07). We will show that it overestimates also rjcfi (but less dramatically). 
The squared matrix element in (|14p for free ultra-relativistic particles can be written as 

(l^-P) = ,2 .2 2 ^^ (18) 



'C 



(^-7), (19) 



where and 7 are dimensionless functions. 



(mfc^-feV/4)(m:V-^V/4) 

^" = WT?)^ ' 

(pL - ^'gV4)(p|, - Y?q^l^) COS^ ^ + h^{v\^ -f pIM^I^ 2 

= n\q^ + A6) ' (22) 

O v/(pk-ft^gV4)(p^,-;iV/4) , ,24 , 

^^11 - -2 ,4(^2 + ^2)(^a+Aa) -c-.c , cos0, (23) 

(j> being the angle between the vectors pi + py and P2 + P2' ■ The function 7 describes interference between two 
scattering channels of identical particles. In the weak-screening approximation, its contribution is small; see Sec. Ill Cl 
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B. Effective collision frequencies 



The collision frequencies are obtained by calculating the integrals ([9]) and pO|) . The calculations are greatly 
simplified because all particles are strongly degenerate. It is sufficient to place the colliding particles on their Fermi 
surfaces (whenever possible) and use the standard energy- angular decomposition based on d'^p = m* pp de dfl, where 
dfi is the solid angle element in the direction of p. All (but one) energy integrations can be done with the aid of the 
energy-conserving delta- function in (|14p : only the lu integration is left. Three angular integrations out of eight are 
performed with the aid of the momentum-conserving delta-function; three integrations (over the position of pi and 
over the azimuthal angle of P2 with respect to pi) are trivial and give Stt^. As a result, one can reduce the angular 
integration to the integration over dq and dcj). Then the collision frequencies ([9]) and (fTO|) can be written as 



w'^ exp(— w) 



7r2p|,^m*c2^ Jo [i_exp(-u.)] 

,,2 



{kBTY 



dw- -I^^i{uj), 



[1 — exp(— w)]^ 



(24) 
(25) 



where w = hco/ (ksT). The functions /f2ci(<^) and I^di^) are the angular integrals 



Inci — 



dq 
"'0 

9 



h^q^ 



2 I 2 
9 [PFc 







dq d0 q^ sj (pIc - f^^Q^ /^){PFi ' ^^974) cos0 ip, 



(26) 
(27) 



where hqm = min{2pc7 2pi} is the maximum momentum transfer in a collision event. Owing to a trivial integration 
over 0, I(ici contains two terms coming from and (p_i_, while J^ci contains only the contribution from f±\\, 



Icic 



(28) 



Let us calculate the angular integrals in the leading approximation with respect to the parameters A/^m and qi/qm- 
This approximation is always justified for the transverse interactions because of the presence of a small quantity 
hu! ~ ksT in the expression for A/q^- However, it is less accurate for the longitudinal contribution since qi/qm is not 
too small; we will discuss corresponding corrections in Sec. Ill CI 

The leading-order expressions for the angular integrals are 



J2 



4 2 
PPcPPi^ 



ml m* pp^c 



2qih^ 



(29) 
(30) 
(31) 



Note that the leading-order expression for is independent of w, being of the same order of magnitude with respect 



to qi/qm as I^^ci- contrast, /j 



r^l'^ (X A" 



The final integration over w gives the collision frequencies 



pV 



he 



1/3 



N5/3 



fl^P^c^l 

2T:^a^'m*pp^ 

fl'^PpcIl 



{kBTf 



(32) 
(33) 
(34) 



where ^ = 2r(8/3)C(5/3)(4/7r)i/3 ^ g 93^ (-(^^^ jg 

the Riemann zeta function, and r(z) is the gamma function. 
Equations ([7]), pT|l - (|13p . and give r}c^ in the weak-screening approximation. 

For typical conditions in a neutron star core, 



K, < K 



K, < K 



(35) 
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However, the inequality is not so strong as for the thermal conductivity (SY07). The dominance of over is 
determined by the factor [hcqi/ {ksT)]^^'^ which increases slowly with decreasing T. It is more accurate to include all 
components of the collision frequencies. For the thermal conductivity problem, we had t'^/i^Jj c>c hcqi/ {ksT), so that 
transverse interactions dominated at all temperatures of interest (SY07). 

Nevertheless, for the not too high temperatures (see Sec. llVl for details), r/o^t is mainly determined by the collisions 
via the exchange of transverse plasmons. In this case, the electron and muon momentum transports are decoupled 
[see Eq. (HH)], 



<^incq,f/^ksTf/\ (36) 
Then the shear viscosity of electrons or muons (c = e or ^) becomes 

We see, that in the low-temperature limit, rj^^ has a non-standard temperature behavior, rjc^ oc T^^^^ (instead of the 
standard Fermi- liquid dependence r] oc T^^). The non-standard behavior was pointed out by Hciselberg and Pethick 
[Sj] for an ultra-relativistic quark plasma. Our results involve collisions of charged particles in the npe/^-matter (for 
any degree of relativity of muons) . Our expressions for r/o^ depend only on the number densities of charged particles 
and on their effective masses; therefore, they can be used for any EOS of dense matter. Previous calculations 0, [isj 
overestimated ?7e^ because they employed the improper plasma screening of transverse interactions. Equation (|37p 
remains valid in the presence of other charged particles (such as hyperons). 



C. Corrections to the leading terms 

As mentioned above, the corrections to containing higher-order powers of h/qm can be neglected, so that 
can be taken in the form (j36p . In contrast, the corrections to vi containing higher-order powers of qi/qm can be 
important. At not too small temperatures, at which can give a noticeable contribution, such corrections can affect 
rjc^- We will discuss several corrections of this type. 



1. Kinematical corrections to vl, and z^' 



The main corrections to the leading terms arise from the g-dependence of the functions (p [Eqs. ([221) and and 
from the (^-dependence in Eqs. ((^ and ((7f)) . The integral J^^- is calculated precisely. 



>k9 >^9 4 9 9 



■lUfin-i/qi) 



4^2 



[mfmfc' +pUrnl' + mf)] l!{q„Jqi) 



+ ^ [«' + m + pIc] 4 ilrn/qi) 



nfi^qf 
64 



4 (9™/®), 



where we have introduced the integrals 



„/k 



■dx'. 



(38) 



(39) 



whose expressions are given in Appendix [X] After the energy integration the corrected collision frequency becomes 



(40) 



Similar corrections should be calculated for In the leading-order approximation (with respect to A/g,„), A^ can 
be neglected in the denominator of f±\\- The remaining angular integral is taken. 



h^qi 



PfcPfzIq iQ^n/qi) - -^{PFc+PFi)Qii2 ii^n/qi) + j^qf {qm / qi) 



(41) 
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where 

After the energy integration we finally obtain 

Our calculations show that these corrections to and v'^^ can reach ~70%. It is advisable to include them in ?7c;i- 



^. Corrections to lepton-proton collision frequencies 



So far we have considered the function ip calculated for free relativistic particles. It is a good approximation for 
collisions within the electron-muon subsystem, because the electrons and muons constitute almost ideal Fermi gases. 
However, the protons belong to a strongly interacting Fermi liquid; this case should be analyzed separately. First of 
all we notice that the protons are non-relativistic. Moreover, we will assume, that many-body effects can be treated 
by introducing an effective proton mass w*. Under these assumptions, the proton transition current can be written 
as Jp2'2 oc ^{P2 +P2')"^p ^(72(^2' ' 'which only slightly modifies (p. The expression (|23p for remains the same, while 
the two other functions become 



CP (pk-^V/4)b|.p-/»2gV4)cos2</, + (p|,^-;,2^V4);iV/4 , 

= ¥WTa^) • (^'^ 

The difference between vlp, calculated with (|44|) and ((2T|) . is proportional to some power of hqi/(m*^c). Contrary 
to qi/qm, this ratio is always small for the conditions in neutron star cores. Hence Eq. (1401) with the angular integral 
psp remains a valid approximation. For the completeness of our analysis, we present the modified angular integral, 

+ ^(.<''+Pfj4(<i.../«)- («) 

This expression gives almost the same Vcp- 



3. Interference corrections to Vcc 



The last correction to be discussed is the correction to i>cc (to v^c and f^i^) due to the interference between two 
scattering channels (1,2 l',2' and 1,2 ^ 2', 1') for collisions of identical particles. This interference is described 
by the dimensionless function 7 which (like ip) contains longitudinal, transverse and mixed components. An accurate 
consideration shows that all these components are smaller than corresponding components of ip. Therefore, the 
interference correction to i/^ can be neglected, and noticeable corrections can arise only to fUc and i/^^. These 
corrections have been calculated in the same way as in previous sections. We have obtained that they are numerically 
small (give < 5% contribution to Vcc + v'^c)- Their contribution to 770^1 is always negligible as expected without any 
numerical calculations. Such corrections can be significant under two conditions. First, the longitudinal component 
Vcc should be comparable to z^^. Second, Vcc itself should give a noticeable contribution to yyo/j- The former condition 
would be realized at high temperatures if particles c are weakly-relativistic. The electrons are ultra-relativistic in 
neutron star cores and do not obey the above requirement. The muons can be weakly-relativistic there, but if they 
are their contribution to 770^ is not large. The importance of the interference corrections in cc collisions is further 
reduced by a (typically) stronger contribution of cp collisions. There are also collisions between electrons and muons. 
The interference corrections for such collisions are absent; corresponding partial collision frequencies are of the same 
order of magnitude as Ucc- 
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Thus, the corrections to Vcp and i^cc seem to be neghgible. The kinematical corrections to i^jj,- and i/^j are significant 

if j^Hj cannot be neglected in comparison with i^^. Note that in SY07, for the electron- muon thermal conductivity 
problem, no corrections have been required because the thermal-conduction fi'equencies ly^^ dominate at any density 
and temperature of practical interest. 



4- Comparison with exact solution 



So far we have used a simplest variational solution for the shear viscosity based on the expression for the trial 
function ([5|) with Tc independent of the particle energy Sc- Actually, however, the energy dependence of <i>c is more 
complicated, which affects the shear viscosity. It is convenient to introduce a correction factor C that relates the 
exact and variational shear viscosities. 



^cxact — V 



(47) 



In ordinary Fermi-systems, where the collision probability is independent of energy transfer huj, the factor C can be 
calculated using the theory developed by Sykes and Brooker for one component systems and extended by Flowers 
and Itoh Q and Anderson et al. for multicomponent systems. Unfortunately, this theory cannot be directly 
applied to our case because of the dynamical character of transverse plasma screening (Landau damping). 

The factor C for the thermal conductivity with account for the exchange of transverse plasmons was estimated in 
SY07. Let us do the same for the shear viscosity. As in SY07, we restrict ourselves to the exchange of transverse 
plasmons in the weak-screening approximation. Then the electron transport decouples from the muon one, and we 
can consider one type of momentum carriers. We redefine <I>c in ([5]) as 

$c = -Toff [vcaVcfi - ^SapVcPc] Va/3^'(a;), 



^^^^ 

where Tog is an effective relaxation time (that can be treated as a normalization constant) , and an unknown function 
'^{x) oi X = [Ec — f^c)/ikBT) describes the energy dependence of ^c- 

Substituting (j48p into the linearized kinetic equation, one obtains an integral equation for ^'(a;), 



/(x)(l-/(x)) 



dx' ■ 



exp(x' — a;) — 1 



f{x){l-f{x')) 



X < ^ ^Inciix' - x)[^{x) - ^{x')] + Indx' - x}^{x') 



where Iqc = ^nd, hid = Inciij and 



dq j d(f) if. 

IQ 



(49) 



(50) 



The integral equation (HS)) is more complicated than that for the thermal conductivity (see Eq. (42) in SY07). The 
term with IqcI in (|49p appears because we go beyond the simplest variational approach of = const. Without that 
approach, there is no cancellation of zero-order expansion terms (in series of q) in kinematical factors in Eq. It 
was that cancellation which led to the appearance of the term under the integral in Eq. (j26p . The integral IqcH 
coincides (save constant factor) with the angular integral for the thermal conductivity problem (Eq. (25) in SY07). 
Taking the weak-screening expressions for /oc and Inci with the exchange of transverse plasmons alone, and choosing 



2/3 



it) aqtihcqt)^/^kBT)5/^' 



(51) 



we obtain the dimensionless equation 



1 



1 + exp(— x) 



dx' 



sgn(a;' — x) 



(exp(a;' — x) — 1)(1 + exp(— x')) 

— oo 

X \x{^ix) - *(x')) + \x' - xp/3^(a;') 



(52) 
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where A = p^p^/ {2>h^ K^). The quantity At = {nkBT/{4Hcqt)y^^qt is the transverse screening wavenumber A, with 
Huj replaced by fc^T. In a neutron star core, one typically has HAt <C Pfcj smd hence A 3> 1. If A ~ 1, then the 
weak-screening approximation is not justified. 

Once ^'(x) is found by solving Eq. ([5^ . the shear viscosity is given by 

oo 

r?coxact = ^^^ll^ / - fix)). (53) 

— oo 

We have solved Eq. ([52|l numerically and compared the result with the variational one (5'var — 2^/^7r^/^/[12^]). For 
A = 10 — 1000, we obtain C = T?exact/??var = 1-08 — 1.056. For A ^ 1000, the factor C ~ 1.056 becomes nearly 
independent of A. Therefore, we have C ~ 1 indicating that the simplest variational approach is well justified. 

Were the electron (and muon) collisions determined solely by the exchange of longitudinal plasmons (with the 
transition matrix element independent of ut), one could find C from the standard theory. In that case one also obtains 
C « 1 (see, e.g., Ref. [IBl)- We expect, that in the most general case, when electron and muon collisions are governed 
by the exchange of transverse and longitudinal plasmons, the correction factor C differs from C = 1 by < 10%. If so, 
the simplest variational approach is sufficiently accurate, and no corrections are required (C = 1) for the majority of 
astrophysical applications. 



D. Neutron viscosity 



In this section we calculate the neutron shear viscosity rjn. We employ the same formalism as was used in BHYOl 
for studying the neutron thermal conductivity. Similar approach was used by Baiko and Haensel to determine 
kinetic coefficients mediated by neutron-neutron collisions. We calculate ?7nn from Eq. the effective relaxation 
time of neutrons, t^, is given by Eq. pip , being determined by the collision frequencies (jl2p of neutrons with neutrons 
and protons. 

The neutron-neutron collision frequency can be written as 

i^nn + ^nn- ^^2^^, ^nn (54) 

(note that the authors of BHYOl did not separate t'nn and v'^^ but considered their sum). The neutron-proton collision 
frequency is 

_ 32m;^ml{kBT)^ ^ 

Here, Wn is a bare nucleon mass. The quantities Snn and S'np are the effective nucleon-nucleon scattering cross-sections 
introduced in Eq. (22) of BHYOl (for the thermal conduction problem). For the shear viscosity, in the same notations 
as in BHYOl, we obtain 



Snn = -rir^ / d.T' / do: , , Q„„, 56 



^2 r0.5+x„ i-a p_(^2 



IGn^h'^ Jo.,.,„ Jo ^ 



Snp = TTrkA / / Ax ' Qnp, (57) 



where x = hq/{2ppn), x' — hq' /{2ppn), a = \/x^ — (0.25 + Xq — x'^)^/^', and xq — pFp/{2pFn)- This choice of 
integration variables is convenient for numerical integration. The quantities Q„n and Qnp are squared matrix elements 
for nucleon-nucleon scattering (in the notations of BHYOl, Qnn = (|AfnnP) and Qnp = (|Mnpp)). 

Let us emphasize that kinematic restrictions in Eqs. ([56|) and ([57|) are very different. The effective cross section 5nn 
is determined by a wide spectrum of momentum transfers q (or, equivalently, of scattering angles). Our calculations 
show that one can get a reasonably accurate S^-a assuming that Qnn is independent of q. Such a g-averaged Qnn can 
be extracted from a total neutron-neutron scattering cross section. In contrast (because, typically, pp-p PFn), Snp 
is determined by small momentum transfers, that is by a small-angle cross section of neutron-proton scattering. 

By analogy with BHYOl, we can write Snn = S^Knn and S'np = Snp Knp- Here, Snn and S'np are the effective cross 
sections calculated with in-medium squared matrix elements, Qnn and Qnp; S^ and S'np^ are similar cross sections 
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calculated with the in- vacuum matrix elements; K^n and K^p are the ratios of the in-medium to the in- vacuum cross 
sections. 

The authors of BHYOl calculated all these quantities for the thermal conduction problem. The squared matrix 
elements Qnn and Qnp were extracted from nuclcon-nucleon differential scattering cross sections calculated in Refs. 
fvA [l^] for symmetric nuclear matter with the Bonn nucleon-nuclcon interaction potential using the Dirac-Brueckner 
approach. An accurate extraction of the in-medium Qnn and Qnp required the knowledge of effective masses m* and 
m* (not reported in T?, Ts'l). For that reason, the procedure used in BHYOl was ambiguous. Thus, the factors K^n 
and Knp, presented in BHYOl for the neutron thermal conductivity, are model dependent and not very certain. 

Now we turn to calculating Snn, Snp, Snn , Snp , Knm and K^p for the shear viscosity. To avoid the above drawbacks, 

we suggest to neglect the in-medium effects on the squared matrix elements and set K^^ = K^p = 1, i^nn = Snn , and 

Snp — . According to BHYOl (for the thermal conductivity), K^n and if^p are indeed ^ 1 (and K^p is relatively 
unimportant). Recently Zhang et al. have studied nucleon-nucleon scattering cross sections in nuclear matter 
taking into account two-nucleon and three-nucleon interactions. They used the Brueckner-Hartree-Fock approach 
and the Argonne V14 nucleon-nucleon interaction model supplemented by three-nucleon interactions. Their principal 
conclusion is that the in-medium effects on square matrix elements are relatively weak, while the main medium effect 
consists in modifying (mostly reducing) m* and m*. The reduction of effective masses under the simultaneous effects 
of two-nucleon and three-nucleon forces is much stronger than under the effect of two-nucleon forces alone. 

Thus, we have calculated Snn and Snp from Eqs. ([56|) and ([57|) using the in- vacuum matrix elements Qnn and Qnp 
from Refs. [13, [3. These matrix elements accurately reproduce [19] well elaborated laboratory measurements of 
differential nucleon-nucleon scattering cross sections. Our calculations of Snn and iSnp are expected to be very close to 
those done with in-vacuum cross sections measured in laboratory. In this sense, our values of Snn and Snp are model 
independent. Similar calculations of Snn in Ref. fld\ give slightly different results due to the different data sets for 
Qnn- Because our equations (|54p and (|55p for the nucleon-nucleon collision frequencies contain a proper dependence 
on m* and to*, they can be regarded as independent of any specific model for nucleon-nucleon interaction. It can be 
a two-body or two-body plus three-body interaction; its explicit form is not essential. Thus, we obtain a description 
of the neutron shear viscosity valid for any EOS of nucleon matter. This approach is not strict (uses the in-vacuum 
matrix elements) but universal. One can in principle calculate more accurate in-medium matrix elements for any 
chosen EOS but loosing the universality. For calculating the diffusive thermal conductivity from the equations of 
BHYOl, we would recommend to adopt the same approach and set i^nn = Knp = 1 in those equations. 

The results of our calculations can be fitted by the expressions 

12.88 1 - 0.6253fcn + 0.3305fc2 

^nn - ^^915 1 - 0.0736fcn ' 

0-8876 fc3-5 1 + 139.6fcp + 103.7fcn 

~ kl 1 - 0.5932fcn + 0.1829fc2 + 7.629fc2 - 0.5405/cpA:n ^ ' ^ ' 

where ki is the Fermi wavenumber of nucleons i expressed in fm^^. As in BHYOl, the calculations and fits cover the 
range of kn from 1.1 to 2.6 fm~^ and the range of kp from 0.3 to 1.2 fm~^. These parameter ranges are appropriate 
to neutron star cores at 0.5 Po ^ P ^ 3po- The fit errors for ^nn do not exceed 0.5%. The maximum fit error of Snp 
is (5max ~ 8% (at kn = 1.1 fm"i and kp = 0.7 fm"^). 



III. SHEAR VISCOSITY IN SUPERFLUID MATTER 

Neutrons and protons in the npe/i-matter of neutron star cores can be in superfluid state (e.g., Ref. Q). Here, 
we study the effects of superfluidity on the shear viscosity. Let Tcnip) be the critical temperature for superfluidity 
of neutrons, and Tcp{p) be the same for protons. Proton superfluidity means superconductivity. Calculations of 
superfluid critical temperatures are complicated and very sensitive to a chosen model of nucleon-nucleon interaction 
and a method to employ many-body (polarization) effects Q. Numerous calculations give drastically different Tcn{p) 
and Tcp{p). It is instructive not to rely on any particular model but treat Ten and T^p as free parameters varied within 
reasonable limits (Tc < 10^^° K), in accordance with microscopic calculations. 

Neutron superfluidity has no direct effect on the shear viscosity yyc^ of electrons and muons (because rjcf^ is limited by 
electromagnetic interactions). However, it strongly affects neutron star hydrodynamics in a complicated way. It makes 
the hydrodynamics essentially multifluid (with several hydrodynamical velocity fields; e.g., j^fli] and references therein); 
it introduces an entrainment effect (which relates motion of neutrons and protons), creates a very specific spectrum 
of elementary medium excitations (phonons) and associated specific energy and momentum transfer mechanisms. All 
these problems go far beyond the scope of our paper. Therefore, we will neglect the effects of neutron superfluidity 
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(will treat neutrons as normal) but consider the effects of proton superfluidity on T]n (assuming the protons to be 
passive scatterers of neutrons, i.e., ignoring momentum transport by protons). We disregard thus hydrodynamical 
effects of neutron and proton superfluids. 

Microscopically, proton superfluidity manifests itself in rearranging proton states (from normal Fermi-liquid quasi- 
particles to Bogoliubov quasiparticles) and in the appearance of a gap A in the proton energy spectrum near the 
Fermi level (e = /i). 



e^fi + sgniOV^+e, (59) 

where ^ = vpip — Pf)', the presented equation is valid at |^| <ti ^i. 

It is generally beheved, that Cooper paring of protons appears in the singlet ^So state (e.g., Ref. 0]). The temper- 
ature dependence of A, calculated in the BCS approximation, can be approximated as |21| 

^ =VT^ri.456-^ + i:I^V (60) 



where r = T/Tcp. 



A. The effects of proton superfluidity on the electron-muon viscosity 

The effects of proton superfluidity on rj^^ are twofold. First, proton superfluidity affects the plasma dielectric 
function and, hence, the screening of electromagnetic interactions. The longitudinal dielectric function is almost 
insensitive to the presence of superfluidity ,22], while the transverse dielectric function modifies collision frequencies 
j'j^. The frequencies I'^j remain almost unchanged because they are independent of the transverse screening in the 
leading order. 

As in SY07, a collision frequency in superfluid matter (to be denoted as i'^^) can be written as 



15 



^iy^,RHy,r), (61) 



where {i =e, ^) stands for a collision frequency in non-superfluid matter, while R^'iy, r) accounts for the superfluid 
effects (which mainly reduce the collision rate); 

r^{pl,+pl^)/plp (62) 

is a slowly varying function determined by plasma composition. We have r = 1 in the absence of muons, and r > 1 
in the presence of muons, with the maximum value of r w 1.26 in the limit of ultra-relativistic muons. The reduction 
factor Rf~{y,r) for the shear viscosity is, however, not the same as for the thermal conductivity (obtained in SY07) 
and will be calculated below. 

The second effect of proton superfluidity consists in an additional (direct; not through plasma screening) reduction 
of the lepton-proton collision frequencies (such as i^cp and I'pp). This direct reduction is exponential; it can be 
described by the reduction factors 



,-LS 
P 

ll-S 



^Tp" = ^cpRp[y.-r), (63) 



^Ip = ^iRliv)- (64) 

The reduction factors and i?p are not the same due to the difference in longitudinal and transverse plasma 
screenings. Following SY07, it is convenient to introduce the reduction factor i?^^ for the total transverse collision 
frequency = J^i + + v^p, 

= i^^^RUy,r), (65) 

RtMr) = [rR^iy,r) + R^{y,r)]/ir + l). (66) 

Below we calculate R^, i?p, i?^ , and i?^ 
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1. Superfluid reduction of collisions in electron-muon subsystem 



Superfluid reduction of lepton-lepton collisions is governed by the transverse polarization function lit. For the 
conditions in neutron star cores, it is sufficient to use Ht in the so-called Pippard limit (huj <^ PFpUppj ^9 ^ PFp and 



^ » l/g, where ^ 
reads 



f^i^Fp/{kBTcp) is the coherence length). In this approximation, the proton contribution to lit 



n 



(p) 



ifp 



(67) 



where q^^ = iap^J (^^tt), and Q is the response function calculated in Ref. [1^ and discussed in SY07 in more details. 

In the non-superfluid limit of ?/ <C 1 one has Q — iTrhw/y, which corresponds to the standard Landau-damping 
expression. In the opposite case of strong superfluidity (y 3> 1), the response function Q becomes pure real, Q = ir^. 
For intermediate superfluidity, y ~ 1, we have used the expressions for Q derived in [23]. They are valid for a pure 
BCS formalism neglecting collective modes and related vortex renormalization in current operators due to gradient 
invariancc. However, as in SY07, the main contribution to 7]^^ comes from the parameter values far from characteristic 
frequencies of collective modes (far from u vpq) and we can use the standard BCS theory. 

The expression for the total polarization function in the superfluid case takes the form 



Ilt = ^{q%^^Q{w,y) 
4qc I- TTW 



2 

9iM 



Qt —'<sQ{w,y) I 

^ TTW J J 



(68) 



In the case of strong superfluidity, the main contribution to lit comes from protons. Moreover, the character of plasma 
screening changes. Instead of the dynamical Landau damping, the screening becomes static, with the frequency- 
independent screening wave number As = [n'^qfpA/{4:Hc)]^^^ . In neutron star cores, one typically has A ~ ksTcp <C 
PFiC- Therefore, the relation Ag <C qi remains true in the superfluid case. In other words, the exchange of transverse 
plasmons in proton superfluid remains more efficient than the exchange of longitudinal plasmons. The strong inequality 
^ Qm justifies the use of the leading-order weak screening approximation in describing the exchange of transverse 
plasmons. 

In order to calculate the reduction factor Rf^ one should reconsider the transverse angular integral I^^^ taking into 
account the changes of electrodynamical plasma properties in superfluid matter. In the leading order with respect to 



where /q ^ refers to non-superfluid matter, while 



'{w,y,r), 



(69) 



F^{w,y,r) 



■Kw{r + l)]^^^ {-Kwr + y'^Q{w, y))^ + {y^Q{w, y)y 



1/3 



X —= sm 
V3 



■ arctan 



\TT'wr + A3(5(w, y)\ 
\-Kwr + y^Q{w, y)\ 



3 y^Q{w,y) 
accounts for superfluid effects. In the limit of strong superfluidity (y 3> 1) we have 



F-^{w,y,r) 



3\/3 



v{r + 1) 
Try 



1/3 



(70) 



(71) 



This asymptotic w-dependence compensates the w-dependence in /^J (that appeared under the effect of plasma 
screening). Moreover, in the expression for /^J the collision energy fko is now replaced by the energy gap A. 
Finally, we write j/^'^ — vjr^Rj^{y,r), and the reduction factor becomes 



(72) 



When superfluidity vanishes (y — > 1) we evidently have R^{y,r) 1. In the opposite case of strong superfluidity 
(y ^ 1) we obtain 

RHy,r) = r^V'^ (73) 

' ' 9V3r(8/3)C(5/3) \ Try J ^ ^ 
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Thus, strong proton superfluidity restores the temperature dependence i/^'^ oc that is standard for Fermi systems. 
This result was derived in SY07 for the thermal conduction problem. It is a natural consequence of changing plasma 
screening from dynamical to statical one when T falls below Tcp. 

In addition, we have computed r) for a wide grid of y. We do not present an appropriate fit, because we will 

calculate and fit the total reduction factor r) for i^^. 



2. Superfluid reduction of collisions of electrons and muons with protons 

Now consider a direct effect of superfluidity on electron-proton and muon-proton collision rates. The consideration 
is similar to that for the thermal conduction problem (SY07). The proton energy gap has to be included in the 
expressions for the collision frequencies through the proton Fermi-Dirac distributions. In addition, the electron- 
proton and muon-proton scattering matrix elements have to be calculated using wave functions of proton Bogoliubov 
quasiparticles. As a result, the reduction factors Rp{y) and Rp{y, r) can be written as 



oo oo 

1 f f dx2 dx2' 



r(8/3)C(5/3) J J l + cxp(z2) 



(Z2' - Z2)|22' - ^21"^/^ (1 + 4U2U2'V2V2') . ,, , , 
7— J TTT ^ F {Z2'-Z2,y,r) 

[1 + exp(-Z2')J[exp(z2' - ^2) - IJ 

(Z2' + -22)12:2' -I- 22!^-^/^ (1 - 4U2U2'V2V2. 



Rliy) 



[1 + exp(z2')][exp(-Z2' - 22) - 1] 

00 00 

3 f f da;2da;2' f (z2' - 2;2)(1 - 4u2U2'V2V2' 



F^i\z2>+Z2ly,r)}, (74) 



r2 







1 -I- exp(2r2) I [1 + cxp(-Z2')][exp(z2' - 22) - 1] 



(Z2' + 22) (1 -I- 4U2U2'V2V2') 

[1 + cxp(z2')][exp(-Z2' - Z2) - 1] 



(75) 



where 



1 / X 



sgn(a;) [ x 

X = VFpip - PFp)/{kBT) and z = {e - fip)/{kBT). 

In the limit of strong superfluidity {y ^ 1) we obtain i?p(y) = cxp(— y) and Rp{y, r) = A-^{r+l)'^/'^y'^/^ exp(— y), 
where 

^ -^io '^'Vo ''^%xp(,?)-exp(,i)"^-^^^2^ 

and 

= —= — / d77i / d?72 -r^^ , w 0.92974. (78) 

3V37ri/3r(8/3)C(5/3) io h exp(772) - exp(77|) ^ > 

Thus, at r <C Tcp collisions with superfluid protons are exponentially suppressed. Then the shear viscosity 770^ is 
limited by collisions within the electron-muon subsystem (which are also affected by proton superfluidity as described 
in Sec, nil A II) . 

We have computed F^iy) for a wide range of y and fitted the results by the expression 
F\(jj) = 1^11 + (1.25- All)cxp(-0.0437 2;) 

-f (1.473 2/2 + 0.00618 ?;^)exp 0.42- ^(0.42)^ + j/2 | 



X exp 



v/(0.22)2+j/2 ^ (79) 
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which reproduces also the asymptotic hmits. The maximum relative fit error is 0.75% at y = 0.533. 

We do not present a separate fit expression for Rp{y, r), but give the fit of the total reduction factor Rj-^^{y, r) 



<t ^ Z ToT? + + 92) exp 0.145 - ^(0.145)2 + ^2 ^ (go) 

(1 + 53r) ^ ^ 

= 0.87-0.314r, 32 = (0.423 + 0.003r)yi/3 + O.OWGy^ - 0.5982/^/^ exp(-y), 
53 = 251r-9(r + 1)6(1 -51)', 

with the maximum fit error ~ 0.3% at r = 1 and y — 3.5. This fit reproduces also the limiting case of R^^t — > 1 at 
?/ — > 0, and the asymptote at y ^ 1, 

47r2r 1 
RtAv,r) = 9^r(8/3)C(5/3)(r+ 1)2/3 ^^^^ 

Recently the electron shear viscosity in superfluid matter has been analyzed by Andersson et al. These authors 
have used the standard (but approximate) approach in which the transverse plasma screening is assumed to be the 
same as the longitudinal one. This approach is inaccurate even in the non-superfluid case. In Ref. 0] the effects 
of superfluidity are described by a reduction factor R*^ for the effective electron-proton collision frequency. That 
factor has been taken from Ref. [l^j devoted to the thermal conductivity problem. However, the reduction factors 
for the thermal conductivity and shear viscosity are different. Moreover, the factor R*^ in Q is inaccurate even for 
the thermal conductivity, because it assumes approximate plasma screening and neglects additional terms associated 
with creation/annihilation of proton Bogoliubov quasiparticles; see SY07 for details. Nevertheless, numerical values 
of 77c, derived from the results of Q for superfluid matter, are not too different from our results. Typically, they 
overestimate rjc by a factor of three, and this overestimation increases with decreasing T^p- 

In the limit of strong superfluidity, the temperature dependence of v^^ formally restores the standard Fermi-liquid 

behavior, vjr^ cx T^. Therefore, vj.^ can be comparable to . Note, that the ratio v'li/i'jr^ in superfluid matter 

remains approximately the same as its value at T = Tcp. If Tcp is sufficiently small, then at T = Tcp we have 3> v^^, 
and the same inequality holds at smaller T . The shear viscosity in superfluid matter, fully determined by the exchange 
of transverse plasmons, can be written as 

Vc = /, — rz^-iz — : ' = jt: -1.71. (82) 



(ksT)^ aV^p^^+pj,^ \pFpcJ ' 40 

One can use this expression for estimates, but we recommend to employ the total collision frequency in practical 
calculations. 



B. Neutron shear viscosity in superfluid matter 



We study the effect of proton superffuidity on neutron-proton collisions. Even this problem is difficult and we adopt 
a simplified approach used in BHYOl for the problem of neutron thermal conductivity. It has also been widely used 
for analyzing superfluid suppression of various neutrino processes (e.g., [25J and references therein). It consists in 
taking an ordinary differential probability of a given scattering process (neutron-proton scattering, in our case) and 
inserting particle energies (j59p with energy gaps in corresponding Fermi-Dirac distribution functions. In our case, this 
approach is expected to be sufficiently accurate. Let us recall that we consider the protons only as neutron scatterers. 
Proton superfluidity suppresses this scattering channel, and our approach reproduces such a suppression. 

In this approximation, the neutron-proton collision frequency becomes 

•^np = I^npRnpiy), (83) 

where i?np(j/) is the superfluid reduction factor. The latter factor is given by the same expression as the reduction 
factor for lepton-proton collisions, i?p, save the coherence factors, 

oo oo 

-P^yi - „2 / / 1 + exp(^2) I [1 + exp(-Z2')][exp(z2' - Z2) - 1] 

(84) 





Z2' + Z2 



[1 -I- exp(z2')][exp(-z2' - ^2) - 1] 
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It obeys the asymptotes i?np(0) = 1 and Rnpiu) ^np2/exp(— y) at y ^ oo, where = 0.8589. 
We have calculated i?„p (y) in a wide range of y and fitted the results by the expression 



Rnpiy) 



0.513+ \/(0-487)2 + 0.018 y2 exp 2.26-^(2.26)2 



1 + 0.00056 y'') exp 



6.2- V(6.2)2 + 4y2 



(85) 



the formal maximum fit error is « 0.25% at y = 11.6. 

Note, that proton superfluidity affects ryn weaker than r]e^ (because of a relatively small contribution of neutron- 
proton collisions to r/n). 



IV. RESULTS AND DISCUSSION 



A. Equations of state 



Our results can be used for a wide range of EOSs of the npe/i-mattcr in neutron star cores. For illustration, we have 
selected five model EOSs. The parameters of these EOSs are given in Table |T1 including the maximum gravitational 
mass Afniax of stable stars and the threshold density of muon appearance. 

Table I; Parameters of the selected EOSs: The compression modulus Ko of symmetric saturated nuclear matter; the muon 
threshold density p^; and also the central density pmax, the mass Mmax and radius -Rm of maximum-mass models (p^ and pmax 
are given in units of lO^** g cm~^) 



EOS 


Ko 


Pm14 


pmax 14 








MeV 






Mq 


km 


APR 


237 


2.28 


27.6 


1.923 


10.31 




120 


2.55 


38.6 


1.468 


9.18 


PAL I 


180 


2.55 


31.4 


1.738 


9.92 




240 


2.55 


26.6 


1.950 


10.59 




120 


2.58 


35.3 


1.484 


9.72 


PAL II 


180 


2.58 


29.5 


1.753 


10.36 




240 


2.58 


25.3 


1.966 


10.97 




120 


2.46 


44.4 


1.416 


8.45 


PAL III 


180 


2.46 


34.5 


1.713 


9.60 




240 


2.46 


28.6 


1.910 


10.12 




120 


2.50 


42.0 


1.438 


8.75 


PAL IV 


180 


2.50 


33.2 


1.713 


9.60 




240 


2.50 


27.8 


1.927 


10.32 



The APR EOS was constructed by Akmal, Pandharipande, and Ravenhall [l^ (their model Argonne 
V18+(5t;+UIX*); it is often used in the literature. Specifically, we adopt its convenient parametrization proposed 
by Heiselberg and Hjorth- Jensen and described as APR I by Gusakov et al. [11]. It is sufficiently stiff, the maxi- 
mum neutron star mass is Mmax ~ 1.92 Mq (and the maximum-mass star has circumferential radius of Rj^ = 10.31 
km), the muons appear at p^ ~ 2.28 x 10^** g cm~'^; see TablelH 

The PAL EOSs are convenient semi-analytical phenomenological EOSs proposed by Prakash, Ainsworth, and Lat- 
timer [29] . They differ by the functional form of the dependence of the symmetry energy S of dense matter on the 
baryon number density nt,. This dependence is described [Hj by a function F{u), where u = Ub/no, no = 0.16 fm~'^ 
being the baryon number density of saturated symmetric matter. For the PAL EOSs I, II, and III, these functions 
are F{u) — u, 2v? /{u + 1), and -^/u, respectively. The PAL IV EOS belongs to the same family of EOSs, but with the 
symmetry energy S{u) oc suggested by Page and Applegate js^]- The PAL EOSs differ also |2§] by the value of 
the compression modulus of saturated symmetric matter, = 120, 180, and 240 MeV. Nevertheless, the particle 
fractions riijn}, as a function of n{, are independent of (for these EOSs); the dependence of nt, on p is almost 
identical for the three selected Kq values [at a fixed F{u)\. Hence, the collision frequencies and the shear viscosity are 
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independent of Kq. However, taking different K^^ one obtains very different neutron star models (different mass-radius 
relations and Mmax; see Table |T]). For illustration, we take Ki^ — 240 MeV for all PAL models (unless the contrary is 
indicated). 

Therefore, our selected EOSs correspond to a large variety of neutron star models. 



B. Shear viscosity in non-superfluid matter 




Figure 1: (Color online) Shear viscosity r^cji of electrons and muons versus density pi4 (in units of 10^^ g cm"^) for different 
EOSs (Table U]) at T = 10* K (rrip = 0.8 mn). The thin solid line (APR-i) shows the viscosity 77c;j calculated with account for 
the exchange of longitudinal plasmons alone. 



Figure [T] shows the shear viscosity rjo^ of electrons and muons versus density at T = 10^ K for five EOSs. The given 
temperature is typical for middle-aged (t ^ 10"* — 10^ yr) isolated (cooling) neutron stars without enhanced neutrino 
emission in their cores (e.g., Refs. [25|,[30|)- The proton effective mass is taken to be m* = 0.8 mn. The thick lines give 
ricf_i for the APR and PAL I-IV EOSs, while the thin solid line is for the APR EOS, but it is calculated including the 
contribution from the exchange of longitudinal plasmons alone. One can see, that the inclusion of transverse plasmons 
lowers 77o/i by a factor of three at p > 4 x 10^"' g cm^'^. With the fall of temperature this lowering is stronger. The 
exchange of transverse plasmons has not been included in previous calculations of the shear viscosity in neutron star 
cores, which has resulted in an overestimation of rjc^. The viscosity rjc^ for the PAL II EOS (the dotted line) goes 
significantly higher than other curves due to a larger amount of protons (and, therefore, electrons and muons) for this 
EOS. 

Figure demonstrates the density dependence of the neutron shear viscosity multiplied by squared temperature, 
ryn T^. This combination is temperature independent. The curves are calculated assuming the nucleon effective masses 
m* = TO* = 0.8 TOn. In principle, the effective masses can be taken from microscopic calculations of an EOS; they can 
depend on p, and our expressions for the shear viscosity allow one to incorporate this density dependence. Here we 
assume density independent effective masses by way of illustration. In Fig. [51 for simplicity, we do not present the 
results for the PAL IV EOS; they are very close to the APR results. One can see, that the neutron viscosity for the 
selected EOSs differs within a factor of < 2. 

Figure [3] demonstrates the viscosity approximation of Cutler and Lindblom [3| (curve CL) versus p at T = lO"^ K. 
Recall that the approximation is based on the calculations by Flowers and Itoh [6i| performed for the EOS of Baym, 
Bethe and Pethick [3l| assuming in- vacuum nucleon- nucleon scattering cross sections and m* = m* = TOu. Also, 
we show self-consistent calculations of ryn by Benhar and Valli [l3| for a pure neutron matter with the EOS that is 
basically similar to APR (with another version for three- nucleon interaction). The authors have used one and the 
same nucleon interaction potential to derive the EOS and ?7„. The curve BV (2bf) is their result (from their Fig. 1) 
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Figure 2: (Color online) Viscosity rjn of neutrons times versus density for four EOSs (m,* — rrip — O.Smn). 



obtained employing two-body nuclcon forces; the curve BV (2bf+3bf) is obtained employing the two-body and three- 
body forces. The effective mass to* is calculated self-consistently as a function of p (m* is different for both curves 
and, unfortunately, is not reported in |1Q|]). 

All other curves in Fig. [3] are our results for the APR EOS assuming various values of to* and m*. For simplicity, 
these phenomenological values are taken density independent. We show either the viscosity rj^^, limited by neutron- 
neutron collisions alone (dashed lines), or the viscosity rji^, limited by neutron- neutron and neutron-proton collisions 
(solid lines). One can see that the contribution of neutron-proton collisions is relatively small, while the dependence of 
the viscosity on nucleon effective masses is important. Smaller effective masses strongly increase the neutron viscosity. 
In the limit of to* = to* = rUn we obtain the viscosity rj^ which is a factor of ~ 40 smaller than CL. Using the results 
of BHYOl for the thermal conductivity of neutrons derived in the same approximations as our results for 77„, 
we obtain the values of a factor of 2-4 smaller than those given by Flowers and Itoh [s^l . The nature of this 
systematic disagreement of our results with the results of Flowers and Itoh is unclear. We have checked that it cannot 
be attributed to using different EOSs. 

A comparison of our results with those of Benhar and Valli [10] is complicated because Benhar and Valli do not 
present the values of to* which they obtained for a neutron matter. If, however, we take a reasonable value of 
TO* = 0.7TOn, we obtain r]nn (not shown in Fig. [3|) very close to the curve BV (2bf) of Benhar and Valli. In order to 
reproduce their BV (2bf+3bf) curve with our equations, we should employ a density dependent to*. It should vary 
from m* « 0.6 mn aX p ^ 1.5 x 10^"* g cm^'^ to 0.45 TOn at p ~ 6 x 10^** g cm^'^. Let us note, that the inclusion 
of three-nucleon interactions does reduce to* and the reduction increases with density [l3|. However, since we do 
not know exact values of to*, used in Ref. |10j . we cannot analyze the relative importance of to* and in-medium 
corrections to the squared matrix element. We assume (Sec. Ill Dl) that the effect of the effective masses is more 
important. Notice, in addition, that the in-vacuum differential neutron-neutron scattering cross section in Ref. [lol | 
(the solid line in their Fig. 3) seems underestimated. 

In Fig. m we compare partial shear viscosities in a neutron star core with the APR EOS at T = 10* K. Previously, 
it has been widely thought that rj^ completely dominates over T^e^i in the core of a nonsuperfluid star. Now we have 
considerably lowered both viscosities (Figs. [T] and [3|). The main contribution to the total shear viscosity (yytot, the solid 
line) at T = 10* K comes from the electrons (?7c, the dashed line). The neutron viscosity rja (the dash-dotted line) is 
lower than ?7c. Note, however, that the relation between 77,1 and rjc is temperature-dependent; when T decreases, "q-n 
becomes more important (see Fig.[5]and a discussion below). The dotted line in Fig.[4]shows the muon shear viscosity 
77^. For T = 10* K, it becomes comparable with r^n at p > 7 x 10^*^ g cm~'^. 
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Figure 3: (Color online) Neutron shear viscosity versus p at r = 10^ K. The curve CL is the approximation of Cutler and 
Lindblom 2] of the results 6]. The curves BV (2bf) and BV (2bf+3bf) are obtained by Benhar and Valli 110] for pure neutron 
matter taking into account two-body and two-body plus three-body forces, respectively. Other curves are our results for the 
viscosity r^nn, limited by neutron-neutron collisions alone, or for the viscosity ?7n, limited by neutron- neutron and neutron-proton 
collisions; these curves are calculated for the APR EOS assuming various phenomenological density independent effective masses 
mn and m*. 



C. Shear viscosity in superfluid matter 

Now we discuss the shear viscosity in the presence of proton superfluidity (superconductivity) but for nonsuperfluid 
neutrons. For illustration, we take m* = m* = 0.8 m,! throughout a neutron star core. 

Figure [5] demonstrates the temperature dependence of rjc^ (solid lines) and (dashed lines) in the presence of 
proton superfluid (Tcp = 10^ K, curves SF) and for non-superfluid matter (unmarked curves) at p = 4 x 10^* g cm"'^. 

In a non-superfluid matter, r/o^ exceeds at T > lO'' K (for the adopted values of to* and to*) but the situation 
reverses at lower T. The reversal is a consequence of the different temperature behaviors, rjc^t T^^/'^ [Eq. (|37p ] 
and cx T~'^. The thin solid line shows ric^ calculated taking into account the exchange of longitudinal plasmons 
alone. It demonstrates the standard Fermi-system behavior, rjcfi oc T~^, and overestimates 770^1 • At T = 10^ K the 
overestimation is small. It reaches a factor of ~ three at T = 10* K, and exceeds one order of magnitude at T < 10^ K. 

Proton superfluidity noticeably increases r^e^ at T < Tcp and restores the Fermi-liquid temperature behavior, 
77c^ cx T^^ [Eq. (HH)]. The increase of jyn is not large because it comes from superfluid suppression of neutron-proton 
collisions which give a relatively small contribution to -q^. In the presence of proton superfluidity, rjc^i completely 
dominates over t^h. 

The electron shear viscosity in superfluid matter has recently been considered by Andersson et al. Q. We have 
already discussed their approach in Sec. IIII A~2l For their EOS and superfluidity model, their results overestimate rye, 
typically, by a factor of three in superfluid matter and by more than one order of magnitude in non-superfluid matter. 

Figure [S] compares the shear viscosity 77 (left panel) with the bulk viscosity C, (right panel) determined by the direct 
and modified Urea processes in the core of a vibrating neutron star at T = 10^, 10®, and 10® K. The vibration 
frequency is set to be w = 10^ s~^; these vibrations strongly affect C but do not affect -q. The bulk viscosity is 
calculated according to Refs. [Hil]. The EOS is the same as in [HlMLiPAL I with Kq = 180 MeV). 

Note, that 77 decreases with growing T, while C increases (e.g., Refs^MUl). For T = 10^ K, the shear viscosity 
dominates in the entire stellar core, while for T = 10* K the bulk viscosity C, in the inner core (where the direct Urea 
process is allowed, after the jump of C in the right panel) becomes ~10 times higher than 77. For T — 10® K, the bulk 
viscosity completely dominates in the entire core. The presence of proton superfluidity enhances 77 and suppresses C,. 
The dashed lines in the left and right panels of Fig.lHlshow 77 and C, respectively, in superfluid matter with Tcp — 10® K 
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Figure 4: (Color online) Partial shear viscosities in non-superfluid neutron star cores versus density at T = 10* K (m^ = rrip = 
O.Smn). 



at T — 10 K. Superfluidity makes the shear viscosity more important. Note that shear perturbations in dense matter 
(e.g., associated with differential stellar rotation) are damped by the shear viscosity and can be unaffected by the bulk 
viscosity. Therefore, the shear viscosity can be important for applications even if it is lower than the bulk viscosity. 

Finally, Fig. [7] compares the shear viscosity in the crust and the core of a neutron star. The viscosity in the crust 
is calculated for cold-catalyzed matter T] using the results of Refs. [35|| and (s^. The former paper is devoted to the 
viscosity mediated by electron-ion collisions, while Ref. [s^] deals with the contribution of electron-electron collisions 
taking into account the exchange of transverse plasmons. In the crust, the latter effect is not large. We use the APR 
EOS in the core, and the core is assumed to be non-superfluid. The solid, dashed, and dash-dot lines correspond to 
T = 10^, 10®, and 10^ K, respectively. The viscosity jump at the star crust-core interface (at p — 1.4 x 10^"* g cm^'^) 
is due to the disappearance of atomic nuclei in the core. The nuclei, present in the crust, lower the viscosity owing to 
a very efficient electron-ion scattering. 

V. CONCLUSIONS 

We have calculated the shear viscosity in a neutron star core as a sum of the electron-muon viscosity ryc;^ and 
the neutron viscosity 77^. Calculating the viscosity r^g^, which is mediated by collisions of charged particles, we have 
taken into account the exchange of transverse plasmons (that has not been done before). Our results include also 
the effects of proton superfluidity. They are universal, presented in the form of analytic fit expressions convenient for 
implementing into computer codes for any EOS of nucleon matter in neutron star cores. 

Our main conclusions are: 

1. The exchange of transverse plasmons strongly reduces rjc^ for all temperatures and densities of interest in a 
non-superfluid core. A low temperatures, we have rjcfi oc T^^^^. 

2. The viscosity r^o^ generally dominates over ryn, although ?]„ can exceed ryo^ at T < 10^ K and p < 4 x 10"'^^ g cm~^ 
(for m* ~ m* « O.Smn). 

3. The viscosity rj^ strongly depends on the nucleon effective masses. Typically, it is more than one order of 
magnitude lower, than that calculated in Ref. @ and parametrized in Ref. ^2]. 

4. Strong proton superfluidity significantly increases r^o/j and restores its Fermi-liquid temperature dependence, 
Tjcfi oc T^^. In this regime, 770^ exceeds ij^. 
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Figure 5: (Color online) Electron-muon and neutron shear viscosities versus temperature in the non-superfluid neutron star 
core and in the presence of proton superfluidity (Tcp — 10^ K) for the APR EOS at p = 4 x 10^'' g cm~^ (m,* — nip = 0.8mn). 
Curves SF correspond to the superfluid case, while other curves are for normal matter. The thin solid curve is r^cn calculated 
including the exchange of longitudinal plasmons alone. 
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Figure 6: (Color online) Density dependence of the shear viscosity rj (left) and the bulk viscosity (right) at different temper- 
atures (the values of logj^p T are given near the curves) in a neutron star core with the PAL I EOS (see text). Solid curves are 
for normal matter (N). Dashed curves are for T — 10* K and proton superfluidity with Tcp — lO" K. The viscosity ^ is plotted 
for a neutron star vibrating at a frequency of a; = 10'* s~^. 



5. The shear viscosity r] is comparable with the bulk viscosity C at T ^ 10^ K (for a star vibrating at a frequency 
uj ^ s~^) and dominates at lower T. Superfluidity increases the importance of r/ in comparison with C. 

Our results can be used in simulations of neutron star hydrodynamics, in particular, to analyze the damping of 
internal differential rotation, stellar oscillations, gravitational wave driving instabilities. 

Our results can be improved further in many respects. It would be most important to study the shear viscosity 
problem in the presence of neutron and proton superfluidity in the frame of multifluid hydrodynamics as discussed in 

Sec. nni 
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Figure 7: (Color online) Density profiles of the shear viscosity through a non-superfiuid neutron star (through the outer crust, 
inner crust and the core) for three values of T (log^g 2^ [K] =7, 8, and 9 - solid, dashed, and dash-dotted lines, respectively). 
The left vertical dotted line shows the neutron drip density, the right line is the crust-core interface. 



Nevertheless, even our restricted standard one-fluid formulation is incomplete. Our calculations of ry„ can be 
improved by taking into account the medium effects on the matrix elements of nucleon-nucleon scattering. However, 
we rely on the results of Ref. [l^ that these medium effects are weaker than the effects of nucleon effective masses 
(which we include explicitly). An account for the medium effects on the matrix elements would complicate the 
expressions for rjn (making them non- universal) . 

We have also neglected the effects of strong magnetic field which can modify the shear viscosity. For not too high 
magnetic fields, B < 10^"^ G, which do not affect the plasma polarization functions (e.g., Ref. 37]), the generalization 
of the present results to the magnetic case is straightforward. For stronger fields, the polarization tensor becomes 
anisotropic and the viscosity problem is very complicated. 

The present results are in line with our studies of kinetic properties of relativistic plasma taking into account the 
exchange of transverse plasmons. These effects were studied by Heiselberg and Pethick for ultra-relativistic quark 
plasma. They should be included in all calculations of kinetic properties of relativistic plasmas, particularly in neutron 
stars. For the neutron star crust, the effect was studied in [3a] (thermal conductivity) and [3y] (shear viscosity). For 
neutron star cores, it was analyzed in [9^ (thermal conductivity), [37| (electrical conductivity), and [39*] (neutrino pair 
emission in electron-electron collisions). 
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Appendix A: EXPLICIT EXPRESSIONS FOR ANGULAR INTEGRALS 



Here we present explicit expressions of the angular integrals ll{x), defined by Eq. (|39p . for different values of k, 
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